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Bosonization approach to the edge reconstruction of two dimentional electron systems 

in a quantum dot 

Itaru YanagQ and Susumu Kurihara 

Department of Physics, Waseda University, Okubo, Shinjuku, Tokyo 169-8555, Japan 

(Dated: February 2, 2008) 

We consider the edge reconstruction of electrons in a two dimensional harmonic trap under a strong 
magnetic field. In this system the edge reconstruction occurs as a result of competition between 
electron-electron interaction and confining potential. To describe it, we develop a bosonization 
scheme for two dimensional electron systems. With this method we obtain the excitation spectrum 
and demonstrate that the edge reconstruction occurs when the value of the magnetic field reaches 
a critical value. We also show that the edge reconstruction depends on the number of electrons. 
Additionally, we calculate the third order terms of bosons in Hamiltonian and examine the effect of 
those terms with a perturbation theory. 

PACS numbers: 74.20.Mn 73.63.Fg 61.46.+W 05.10.Cc 



I. INTRODUCTION 

In recent years technological advances have made it 
possible to construct ideal two dimensional electron sys- 
tems in quantum dots Q . These devices allow us to study 
various behaviors of correlated electron systems. 

In this paper, we study the edge reconstruction of two 
dimensional electron systems in a quantum dot under a 
strong magnetic field. It is a transition caused by inter- 
play between electron-electron interaction and confine- 
ment potential. At relatively strong magnetic fields, all 
electrons are spin polarized and exist only in the lowest 
Landau level (LLL). When the magnetic field is below a 
critical value, electrons in a dot form a compact droplet 
(filling factor v = 1), occupying orbitals with the lowest 
total angular momentum. If the magnetic field reaches 
the critical value, the compact droplet is replaced by the 
system composed of a smaller droplet surrounded by a 
ring, with vacancies introduced not into the center, but 
rather near the droplet edge. Such a transition has been 
studied with numerical analyses 0, fj . Our purpose is to 
describe it in more analytical form, using bosonization 
approach. 

There are some other works about bosonization for two 
dimensional electron systems. Westfahl et al. introduced 
Landau level bosonization of two dimentional electron 
systems Q. Rojt et al. reported Tomonaga-Luttinger 
(TL) liquid behavior in a two dimensional circular quan- 
tum dot at zero magnetic field 0- 

The edge reconstruction is also studied at fractional 
quantum Hall edges (FQHE) [1 E3, EH El . The differ- 
ence between FQHE and our model is a shape of poten- 
tial. While the confinement potential of our model is a 
weak hermonic trap, FQHE is made by sharp boundaries 
and supposed to form chiral Luttinger liquids The 
mechanism of the reconstruction is essentially the same, 
i.e., interplay between electron-electron interaction and 



confinement potential. 



II. MODEL AND HAMILTONIAN 

In our model, electrons are confined in a two dimen- 
sional harmonic trap under a strong magnetic field per- 
pendicular to the plane, 
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where TV is the number of electrons in the quantum dot, 
to* is the effective mass of electrons (to* = 0.067m e ), uq 
is a harmonic oscillator frequency and V(|r|) is Coulomb 
interaction of electrons. 

Now all electrons are assumed to be spin-polarized and 
exist only in the LLL due to a strong magnetic field, the 
noninteracting part of the Hamiltonian is easily diago- 
nalized and its eigenstates and energies are given by 




where n is a quantum number of an angular momen- 
tum (n = 0,1,2...), lo c is a cyclotron frequency, Iq — 
\Jhj2m*Q, and £1 = ^/uj^+ujJJZ. The degeneracy in the 
LLL is lifted by a parabolic potential. 



III. COMPARISON WITH ONE DIMENSIONAL 
ELECTRON SYSTEMS 
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In order to compare this model with one dimensional 
TL model, we rewrite the Hamiltonian in terms of second 



2 



quantization, 



j,l c l-k c j+k c 3 c l > 



(4) 



,k,l 



where cj is an operator that creates an electron with 
angular momentum / in the LLL and Vkj i is the matrix 
element of Coulomb interaction. 

e; has the linear dispersion similar to the kinetic energy 
part in TL model. Our model, however, has only one 
linear dispersion branch while TL model has two linear 
dispersion branches, so we don't have to consider the 
interaction between two different branches. 

Vk,j,i has the essential deviation from the TL model 
because it has not only fc-dependence but also j, l- 
dependence. In the case of one dimensional electron sys- 
tems, there is only fc-dependence which corresponds to 
the momentum transfer between two particles. 

If fc-dependence is much stronger than j, /-dependence, 
it would appear that a bosonization method is effective 
for our model. Vkjj can be calculated exactly in this 
model 0, HJ • The Appendix shows the function form of 
Vi,j,k- We evaluate fc-dependence and j, /-dependence of 
Vfcj,;, then we find that fc-dependence is more dominant 
than j, /-dependence ( Vk+ij,i —Vkj,i » Vkj+i,i—Vkj,i, 
less than 1/10 for several values of i, j, fc ). 

But then, the edge reconstruction is caused by j, l- 
dependence, so we have to develop a bosonization scheme 
including j, /-dependence. 



IV. BOSONIZATION FOR TWO 
DEIMENSIONAL ELECTRON SYSTEMS 

First, we describe fermionic operators c, in terms of 
bosonic operators b, b\ Bosonization of fermionic oper- 
ators has been studied well Q, but this time we cannot 
use these method without modification. We follow the 
prescription proposed by Schonhammer and Meden [6|. 
Creation and annihilation operators of electrons can be 
described with boson operators as follows, 



terms of boson operators, we can obtain a bosonized form 
of the Hamiltonian, 
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here, the ground state is the vaccum of bosons, fc is 



defined as fc 



where £ m i m is a total angular 



momentum of electrons which occupy the orbitals with 
the lowest values of angular momentum allowed by the 
Pauli principle. Vkjj equals when j or / is negative. 
Higher order terms of bosons are less dominant when 
N >> fc. The coefficients of third order terms of bosons 
are much smaller than those of quadratic order terms 
(less than 1/10 of harmonic terms) as will be described 
in detail later. We consider only quadratic order terms 
of bosons in this section. 

Comparing the Hamiltonian (8) with the bosonized 
form of TL Hamiltonian, the first term in the coefficient 
of b\pk has the same form as that of TL Hamiltonian. 
The second term is also almost the same and exactly cor- 
responds to that of TL hamiltonian when j, /-dependence 
does not exist. The third term is the most important in 
this case because this term causes the edge reconstruc- 
tion. It is always negative, so there is a possibility that 
the excitation energy also becomes negative. This term 
is does not exist without j, /-dependence. 

Additionally, the Hamiltonian (8) has several remark- 
able features. It includes the terms of Voj : i (angular mo- 
mentum transfer=0) and N which are ordinary classified 
into the ground state in TL model. 

Fig.l shows the excitation spectra for several values of 
magnetic field B when N = 50 and Huiq = 2.3[meV]. We 
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where ifi(v) is an auxiliary field, N is the number of elec- 
trons and b n is an operator which annihilates a boson 
with an angular momentum n. Substituting these equa- 
tions into eq.(4) and expanding it up to quadratic order 



I 

c 
<a 

p 

5 



B=5[T] t 



****** B=6[T] 
************ 



! ! ! 



■ . b B=S.5[T] 



B 

2 4 6 8 10 12 
angular momentum k 

FIG. 1: Excitation energy spectra for several values of the 
magnetic field, huio = 2.3[meV] and = 50 
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FIG. 2: Excitation energy spectra for several values of the 
number of electrons, huo = 2.3[meV] and B = 8.5[T] 
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FIG. 3: Occupation coefficients of single particle orbitals 
< c„c n > considering up to quadratic order terms of bosons 
in eq.(8). This graph shows the occupancies when a boson 
with an angular momentum k = 11 is excited. The range of 
a reconstructed part is from n — N — k + 1 to n — N + k and 
expressed as a simple step function. 



take the number of electrons sufficiently large to be able 
to use the bosonization method. This figure shows that 
the edge reconstruction occurs with increase of magnetic 
field. As it increases, the difference between kinetic ener- 
gies of neighboring single particle orbitals decreases, and 
the energy of the coulomb interaction increases. Those 



are the causes of the edge reconstruction. In our work, 
the edge reconstruction occurs between B — 8.5[T] and 
B = 9.5 [T]. Moreover we can observe a characteristic 
minimum of the energy at k = 13 for B = 6[T], such 
minimum is also confirmed in the numerical diagonaliza- 
tion and referred to as the edge magneto-roton [j] . 

Fig. 2 shows the energy spectra for several values of 
the number of electrons when Hluq — 2.3[meV] and 
B = 8.5 [T], We confirm that the edge reconstruction 
occurs with decrease of the number of electron N, due to 
the increase of the energy of Coulomb interaction. The 
dependence of the energy spectra on the number of the 
electrons is also investigated with the numerical diago- 
nalization 4], and the behaviors of the energy spectra 
in this paper is consistent with those in the numerical 
analysis. 

Next, we calculate the occupation coefficients < 
Cn°n > °f single particle orbitals which can be calculated 
with eq.(5)~(7), 



<k\clc n \k> = 2J(1 



N,-l 



1=0 



+ -j^(Sn-N-l-k + 8 n -N-l + k)] , (9) 

where \k > is the state in which an elementary excita- 
tion of bosons with angular momentum k is excited. We 
cannot get appropriate occupation coefficients as long as 
considering up to quadratic order terms of bosons. Fig. 3 
shows an occupation coefficient and we cannot confirm 
the configuration that a smaller droplet is surrounded by 
a ring, which is one of the important feature of the edge 
reconstruction and confirmed in the numerical analysis. 
So we have to consider higher order terms of bosons to 
obtain such an appropriate configuration. 



V. THIRD ORDER TERMS OF BOSONS 

In this section, we discuss the third order terms of 
bosons in Eq (8). It can be calculated as follows, 
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(10) 



As mentioned previously, most of the coefficients of third quadratic order terms, so we treat third order terms with 
order trems are small enough compared to those of 
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a perturbation theory, we take into account up to second- 
order perturbations with respect to energies and up to 
first-order perturbations with respect to energy eigen- 
states. Both of them include processes that a boson with 
angular momentum fc splits into two particles with angu- 
lar momentum fc — n and n as an intermediate state, 
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Where \k — n,n > is the state that bosons with angular 
momentum k — n and n are excited and \k > pcr t. is the 
state including the effect of third order terms as a pertur- 
bation. [k/2] means the integer part of k/2. E^ is the 
unperturbed part of energy i.e., quadratic order terms 

(2) 

in eq.(8). E k is the second-order perturbation with re- 
spect to Hz r< i- The first-order terms of the perturbation 
does not appear. 

Fig. 4 shows excitation spectra when huj = 2.3[meV], 
B = 9.5[T] and N — 50. The spectrum A is unperturbed 
one and the other spectrum B includes the effect of per- 
turbations. We can confirm the decreases of energies due 
to the second-order perturbation. As an angular momen- 
tum k increases, a deviation from the unperturbed spec- 
trum also increases because the number of the possible 
combinations of k — n and n increases. 

Figs. 5 and 6 show the occupation coefficients < c^Cn > 
calculated with a perturbed state, eq.(12). The function 
form of < c\ x c n > is so complicated, composed of such 
a great number of Kronecker deltas that we don't tran- 
scribe it on this paper. These graphs show the occupan- 
cies of the ground state with total angular momentum 
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FIG. 4: Excitation energy spectra when huio = 2.3[meV], 
B — 9.5[T] and iV = 50. The spectrum A is unperturbed and 
B is perturbed. 
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FIG. 5: Occupation coefficients < cj,c„ > including the effect 
of third order terms of bosons when N = 50 and total angular 
momentum Ltot = 56 (fc = 6). 
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FIG. 6: Occupation coefficients < c„c„ 
of third order terms of bosons when iV : 
momentum Ltot =61 (fc = 11). 



> including the effect 
= 50 and total angular 



Ltot = 61 (fc = 11) and L tot = 56 (fc = 6) when N = 50. 
Occupancies are fixed independent of confinement poten- 
tial and magnetic field because eq.(12) is independent of 
<x>o and B. We can confirm the configurations from both 
graphs that a smaller droplet is surrounded by a ring. 
These configurations shows the tendency that the recon- 
structed part escapes further from the center of the dot 
than in the case of Figure 3. The more vacancy exists be- 
tween the edge and the center, the more Coulomb inter- 
action reduces. This means that the third order terms of 
bosons make a contribution to the reduction of Coulomb 
interaction energy. 



VI. SUMMARY 

In summary, we have developed a bosonization scheme 
to enable to treat such Hamiltonians as (4). Using this 
method, we have obtained the excitation spectra analyti- 
cally for several values of B and N and observed the edge 
reconstruction. We have also observed the appearance of 
the edge magneto-roton. 

In addition to these results, we have examined the ef- 
fect of the third order terms of bosons and obtain ap- 



5 



propriate occupancies. The third order terms of bosons 
produce a ring at the edge of the dot and reduce the 
energy of Coulonb interaction. 

Such behaviors as above agree well with the numeri- 
cal analysis and we can say that the edge reconstruction 
has been recognized in a more analytical form, with a 
bosonization method. 
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There are some other representations of Vk,j,i, e.g., in 
references 0, and which are consistent with ours. 



The matrix element of Coulomb interaction I4,j,fc 
written as follows, 
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